Light propagation in a birefringent plate with topological charge 
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We calculated the Fresnel paraxial propagator in a birefringent plate having topological charge q at 
its center, named "q-plate". We studied the change of the beam transverse profile when it traverses the 
plate. An analytical closed form of the beam profile propagating in the "g-plate" can be found for many 
important specific input beam profiles. We paid particular attention to the plate having a topological unit 
charge and we found that if small losses due to reflection, absorption and scattering are neglected, the 
plate can convert the photon spin into orbital angular momentum with up to 100 % efficiency, provided the 
thickness of the plate is less than the Rayleigh range of the incident beam. © 2009 Optical Society of America 
OCIS codes: 050.1960, 260.1960, 260.6042. 



Light beams carrying orbital angular momentum (0AM) 
are receiving increasing attention as a resource in quan- 
tum and classical optics, since 0AM exists in an inher- 
ently multidimensional space. Information can thus be 
encoded in higher dimensional 0AM- alphabets [1,2] for 
its use in free-space communication systems [3] or to in- 
crease the dimensionality of the working Hilbert space in 
quantum communications systems [4] . The main charac- 
teristics of a light beam carrying 0AM is the presence 
of a topological point-charge of integer order in the op- 
tical phase. It was recently demonstrated that a bire- 
fringent plate made of a vortex-patterned liquid crystal 
film can imprint its topological charge into the optical 
phase of the incident light, thus changing the 0AM of 
the beam. Being birefringent, the plate affects the pho- 
ton spin angular momentum (SAM) too, thus providing 
easy control of the beam 0AM content by either chang- 
ing its polarization or changing the retardation of the 
plate [5]. It can be shown that, if the topological charge 
of the birefringent plate is q, the 0AM of a light beam 
passing through such a "g-plate" (QP), changes by an 
amount ±2qh per photon. One may roughly think of the 
QP as of an ordinary birefringent plate rotated at angle 
a about the beam z-axis, with a given by a = a{x, y) = 
arctan(2//x) = 4> where is the azimuthal angle in the 
X, y-plane. This simple picture is enough to account for 
some qualitative effects such as, for instance, the SAM- 
to-OAM conversion (STOC) [5], the associated optical 
Berry phase [6] and the changes of the SAM and 0AM 
content of the beam at different depths in the QP [7], 
but cannot be used to determine the detailed quanti- 
tative behavior of the beam during its propagation. In 
fact, the formal replacement a — > 4>(x, y) is justified only 
if the function a{x,y) varies smoothly over the optical 
wavelength scale (the so-called geometric optics approx- 
imation (GOA)) [7], which is not true in the present 
case, because no length scale is defined by the function 
arctan(y/a;). As shown below, exact series solution of 



the wave equation inside the QP shows that all field and 
wave propagators at the QP singularity vanish which is 
not true for fields [5] and propagators [7] calculated by 
just setting a ^ </>. To overcome this broad discrepancy, 
in this work, we study the propagation of a light beam 
in the QP without having recourse to the GOA, but as- 
suming, instead, a good beam paraxiality. We start from 
Maxwell's wave equation V^E - V(V-E) + kl e-E = 0, 
where fco = 27r/A ~ lo/c, e is the relative dielectric 
tensor at frequency w, c is the speed of light and A is 
the wavelength. The relative dielectric tensor is given by 
e = R{q(j))- eiocai- Ri-qcj)), where eiocai = <ii&g{nl,nl,nl) 
is the dielectric tensor in the local frame of the plate, 
R{a) is the rotation matrix about the z-axis, </> is the az- 
imuthal angle, Uq and rig are the material ordinary and 
extraordinary indices, respectively. We neglect absorp- 
tion and assume positive birefringence n,, > Uq. In most 
materials, including liquid crystals, the birefringence is 
small (jT-e — Uo) ^ Tio, so that we may neglect the lon- 
gitudinal part of the optical field and take V-E ~ 0. 
In this approximation, the wave equation reduces to 
Hclmholtz's vector equation V^E_l -|- fcge-E_L = for the 
transverse part Ej^ of the field. In view of the cylindri- 
cal symmetry of the problem, it is convenient to find the 
eigenmodes of the Helmholtz's vector equation in the cir- 
cular polarization basis E± = {Ex ±zi?y)/-\/2 and in the 
cylindrical coordinates (r, 4>, z), by setting Ej_(r, 0, z) = 
(£+(r)e*("+9)'^,£_(r)e*('"-'?)'^,0)e-^'=«'^^'+*"*, where 7 is 
the longitudinal spatial frequency and z = is the input- 
face of the QP. Inserting this field into Helmholtz's equa- 
tion, yields a pair of coupled radial equations 

where /(r) = (£+ + £„)/V2, g{r) = {£+ - £-)/V2, 
H = \J -I- q^ and v = 2m q . Equations (1) are exact 



and can be solved by series. We postpone a full discussion 
about their solutions to a future work. Here we consider 
only the approximate solutions for paraxial beams at 
normal incidence. Setting 7 



7o 



/? /2no, where /? is a transverse spatial spectrum, Eqs.(l) 
reduce to the equations for the ordinary wave 
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where = 1 + 2{n1 — n^^)/ 0^ . We observe that in com- 
mercial liquid crystals we have nl — ~ 0.5 while 
usual paraxial laser beams at normal incidence have the 
radial spatial frequency with P ranging from zero to 
f3 ~ 10~^. The parameter is therefore very large in 
all practical cases. We may then solve Eqs.(2) as asymp- 
totic series of A^. The zero-order approximation of the 
asymptotic solution of Eq.(2) for the ordinary wave is 
given by /o(r) = AoJf^{ko(3r) and goir) = 0, where 
Ao is an arbitrary constant and J^{x) is Bessel' function 
of index ji = \J m? + . The differential equation for 
the extraordinary wave is obtained from Eqs.(l) by set- 



ting 7 
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(3 /2ne and it can 



be obtained from Eqs.(2) by the formal replacements 
/ g,g — > /, A^ — *■ — A^. The zero-order asymp- 
totic solution for the extraordinary wave is then given 
by fe{r) ~ and ge{r) ~ Af.J^{kQf3r) with constant 
Ag. All other terms of the asymptotic solution can be 
found recursively. The zero-order asymptotic solutions 
hold in the whole x, y-plane except a small region about 
the origin having a radius tq « X/{Kj3) = \/ ^Jnl ~ n1. 
For commercial liquid crystal we have rg ^ 1.25 A. In 
this small region is located the singularity and the op- 
tical axis is not well defined here. The effect of this re- 
gion can be accounted for only by the exact (not parax- 
ial) wave approach, but we may anticipate on physical 
grounds that the main effect of this region is to scatter 
a small fraction of the light at large angles out of the 
paraxial beam. From the asymptotic paraxial modes of 
Hclmholtz' equation it is straightforward to calculate the 
Frcsncl paraxial propagator for the optical field E_l. The 
optical field at plane z in the QP is given by 
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+ (A-°-A-)a,}i?(-#)Ei(p,7/.,0). (3) 



where 1, R{4>) and az are the 2x2 unit, rotation and 
Pauli's matrices, respectively. The Fresnel kernels in 
Eq.(3) are given by K"^' = ^;(™)(^, P; 2)e™('^"'^\ 
where ^{m) = ■\/m^~4^^ and 
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The Fresnel kernels K° and in Eq.(4) arc character- 
ized by the presence of Bessel function of irrational or- 
der. Although, K" and K'^ cannot be obtained in a closed 
form, they permit to evaluate analytically the field trans- 
mitted by the QP in important cases as, for instance, for 
Laguerrc-Gauss incident beams. Here we consider only 
the case of a LGq; beam impinging onto the QP. Set- 



ting Ei(p,0,O) = e»'^LG[,(p) 



in the circular po- 



larization basis, where LGq(p) is the radial amplitude of 
Laguerrc-Gauss modes, we obtain 
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where = (HyGG|i|_^.^(r, z/rio) ± 

g-.fcoA„.jjyGG|i|_^,^(r,2K))/2, An = ~ n„, 

/i"*" = ^{£±q) and HyGGp ,„(r, 2) is the 
Hypergeometric-Gaussian mode [8], viz. 
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where r - + l / 2p+I"'I + i r(l + |m|-Hf) / 

wtiereOp„-ii y rf(p+|™|+i) r(|mH-i) ^P-^/^o, 
( = z/zn and zn = kQw'^/2 is the beam Rayleigh 
range. Because Uo — rie, the arguments of the func- 
tion HyGGp,„ in Eq.(5) are very close, so that when 
Anz = jX {i = 1,2,...) the matrix in Eq.(5) is al- 
most diagonal, the beam in the QP has the same value 
of 0AM, i.e. tfi per photon. When Anz = {2j - l)A/2, 
instead, only the off-diagonal elements survive, the right 
and left circular components of transmitted field assume 
a phase factor e^^'^*^ and the beam 0AM change by 
±2qh per photon, depending on the input circular polar- 
ization helicity. As the beam propagates in the QP, its 
transverse profile, spin and 0AM change. From Eq.(5), 
we may calculate the average SAM and 0AM carried by 
the beam at the plane z in the QP, obtaining 
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and xiO 



rtpno — tine — n. 



As expected, we have 



/„.„(0) = 1 so that Eqs.(8) and (7) yield 540) = (| ftp 
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Fig. 1. Beating of SAM (dashed line) and 0AM (full 
line) as a function of the optical retardation Anz / A while 
a circularly polarized input beam propagates in the 1- 
platc. (a) For LGoo and (b) LGoi as a input beam. We 
used the following data: Uo = 1-5, Ue = 1.7, wo = 50A. 

\a\^)/LO and L,{0) ^ (jop + \b\^)£/Lu. In Fig.(l) the pho- 
ton STOC [5] is shown as a function of the propagation 
depth in the 1-plate for LGqo and LGqi input beams. The 
conversion efficiency is practically 100% and its maxi- 
mum occurs at optical retardation Anz — {2j — l)A/2 
with integer j. When the optical retardation is jA, no 
conversion occurs and the beam has no 0AM. Changing 
the optical retardation of the 1-plate provides a good 
way to control the STOC process. However, when the 
thickness of the 1-plate becomes very large (much larger 
than the beam Rayleigh range) the conversion efficiency 
slowly decays. According to Eqs.(5) and (6), the field 
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Fig. 2. Intensity profile for (a) fuU STOC (b) no STOC 
in the 1-platc. Full and dashed lines are simulated by [7] 
and our theory, respectively. The input beam assumed 
the TEMoo. 

profile inside the 1-plate (and at its exit face) vanishes 
as along the beam axis so that the intensity profile 
has the characteristic doughnut shape irrespective of the 
0AM carried by the beam. Fig. 2 shows the intensity 
profiles for (a) fuU STOC (b) no STOC. For the sake 
of comparison, the results obtained in the GOA [7] are 
also shown. We can deduce that the GOA approxima- 
tion is fairly good for the case of full STOC, but is very 
bad in the near field and in the case of no STOC. Dra- 
matic changes of the intensity profile depending on the 
final 0AM are seen, however, in the far-field after free- air 
propagation. When the STOC is maximum, in fact, we 
observe the doughnut profile, while when no conversion 
occurs, the far-field pattern has again a maximum at its 
center. This is shown in Fig. (3). 

In conclusion, we calculated the Fresnel propagator of an 
optical beam in a birefringent plate having a topological 



Fig. 3. Intensity profile in the far-field beyond the 1-plate 
after free-air propagation, (a) No STOC; (b) Full STOC. 



charge q in the paraxial approximation and for normal 
incidence. We considered in the some details the propa- 
gation of a LGqi beam in the QP. As the beam traverses 
the plate, its transverse profile changes from Laugerre- 
Gaussian to Hypergeomctric-Gaussian and STOC oc- 
curs. We paid particular attention to the 1-plate and we 
found that the conversion efficiency is almost 100% when 
the thickness of the plate is much smaller than the beam 
Rayleigh range and slowly decreases when the thickness 
of the 1-plate is increased. The free propagation of the 
HyGG modes generated by the QP is not stable and 
the characteristic transition from dot-to doughnut pro- 
file when the 0AM changes from is observed only in the 
far field. The possibility offered by the azimuthally ori- 
ented plate in manipulating entanglement among several 
degrees of freedom of the light may be of great interest 
for quantum information, quantum communications and 
quantum computing. 
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